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Abstract 

We include the full second-order corrections to the static QCD potential in the analysis 
of the ti threshold cross section. There is an unexpectedly large difference between the 
QCD potential improved by the renormalization-group equation in momentum space and 
the potential improved by the renormalization-group equation in coordinate space. This 
difference remains even at a fairly short distance 1/r ~ 100 GeV and its origin can be 
understood within perturbative QCD. We scrutinize the theoretical uncertainties of the 
QCD potential in relation to the ti threshold cross section. In particular there exists a 
theoretical uncertainty which limits our present theoretical accuracy of the ti threshold 
cross section at the peak to be 5<T pe ak/cpeak ^ 6% within perturbative QCD. 



*On leave of absence from Department of Physics, Tohoku University, Sendai 980-77, Japan. 



In this paper we report on our present theoretical understanding of the ti total cross sec- 
tion near the threshold. Up to now, all the 0(a s ) corrections (also leading logarithms) have 
been included in the calculations of various cross sections near threshold. In order to take 
into account the QCD binding effects properly in the cross sections, we have to systematically 
rearrange the perturbative expansion near threshold. Namely, we first resum all the leading 
Coulomb singularities ~ (a s //3) n , take the result as the leading order contribution, and then 
calculate higher order corrections, which are essentially resummations of the terms ~ a™ +1 //3 n , 
a™ +2 / f3 n , ... It is also important to resum large logarithms arising from the large scale difference 
involved in the calculation JT|.[] This is achieved by (first) calculating the Green function of the 
non-relativistic Schrodinger equation with the QCD potential 0, []J. Conventionally both the 
coordinate-space approach developed in Refs. [|T], |3] and the momentum-space approach devel- 
oped in Refs. ||, f| have been used in solving the equation by different groups independently. 
It has recently been found P that there are discrepancies in the results obtained from the two 
approaches reflecting the difference in the construction of the potentials in both spaces. It was 
argued that the differences are formally of 0(a 2 s ) but their size turns out to be non- negligible. 

Quite recently there has been considerable progress in the theoretical calculations of the 
second-order corrections to the cross section at threshold and the Coulombic bound-state prob- 
lem. New contributions have been calculated analytically |7], [8] and numerically for QED 
bound-states. Very important steps have been accomplished in QCD as well. The full second- 
order correction to the static QCD potential was computed in f|3T0| | . Also, the 0(a 2 ) total cross 
section is known now in the region a s -C f3 -C 1 as a series expansion in (3 |PH| . All these 
results have to be included in the calculation of the full 0(a 2 ) corrections to the threshold 
cross section, which has just been completed (as far as the production process of top quarks are 
concerned) |]12[ . The full second-order corrections turned out to be anomalously large, which 
may suggest a poor convergence of the perturbative QCD in the ti threshold region. 

In this paper, we incorporate the full 0{a1) corrections (the second-order corrections to the 
leading contribution) to the static QCD potential into our analyses. In principle this is a step 
towards an improvement of the theoretical precision in our analysis of the ti threshold cross 
section. Then we scrutinize the problem of the difference between the momentum-space and the 
coordinate-space potentials. Contrary to our expectation, the inclusion of the above corrections 
does not reduce the difference of the cross sections significantly, and there still remains a non- 
negligible deviation. We find that there is a theoretical uncertainty within perturbative QCD 
which limits our present-day theoretical accuracy of the threshold cross section. 

Let us first state the numerical accuracies attained throughout our analyses. We confirmed 
that our numerical accuracies are at the level of 10~ 4 . We have tested our programs with the 
Coulomb potential whose analytical form is known both in momentum space and in coordinate 
space. Moreover we confirmed that we obtain the same cross section within the above accuracy, 
irrespective of whether we solve the Schrodinger equation in momentum space or first Fourier 
transform the QCD potential and solve the Schrodinger equation in coordinate space. In this 
way we also checked that our numerical Fourier transformation of the QCD potential (from 
momentum space to coordinate space) works within the above quoted accuracy. The level of 

* Since the toponium resonance wave functions have wide distributions ~ 10-20 GeV, they probe a fairly 
wide range of the QCD potential. For example, this is reflected in the fact that the fixed-order calculation with 
any single choice of scale /i cannot reproduce simultaneously both the distribution and the normalization of the 
differential cross section which includes all the leading logarithms. It is known that the normalization of the 
cross section is more sensitive to the short-distance behavior of the QCD potential. 
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accuracies is quite safe in studying the size of the higher order corrections which are described 
in this paper. 

Let us now briefly explain the construction of our potentials in momentum space and in co- 
ordinate space, respectively. More detailed descriptions including formulas are given in the ap- 
pendices. The large-momentum part of the momentum-space potential Vjkpt(<?) is determined 
as follows. First the potential has been calculated up to 0{a z s ) in a fixed-order calculation. 
The result is then improved using the three-loop renormalization group equation in momentum 
space. At low momentum, the potential is continued smoothly to a Richardson-like poten- 
tial. On the other hand, the short-distance part of the coordinate-space potential Vsfhmn(^) 
is calculated by taking the Fourier transform of the fixed-order perturbative potential in mo- 
mentum space, and then is improved using the three-loop renormalization group equation in 
coordinate space. At long distance, the potential is continued smoothly to a phenomenological 
ansatz. Thus, it is important to note that the two potentials are not the Fourier transforms of 
each other even in the large-momentum or short- distance region. They agree only up to the 
next-to-next-to-leading logarithmic terms of the series expansion in a fixed MS coupling. The 
difference begins with the non-logarithmic term in the three-loop fixed-order correction. 
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Figure 1: Comparison of the total cross sections (normalized to R) calculated from the different 
potentials: Vjkpt (solid), Vsfhmn (dashed), and V new (dotted line). We set a^g(Mf) = 0.118, 
m t = 175 GeV, and r t = 1.421 GeV. 

In Fig. |I] we show a comparison of the total cross sections (normalized to R) calculated from 
Vjkpt (solid) and from Vsfhmn (dashed line), without any weak or hard-gluon corrections: 

4 poo 

R=— dpp 2 \G(E,p)\ 2 T t . (1) 
Tirrit Jo 
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For the physical parameters we used a^g(Mf ) = 0.118, m t = 175 GeV, and T t = 1.421 GeV. 
We find that the two cross sections differ by 7.8% at the peaks and by 2.2% at E = 5 GeV.Q 
Since the difference of the cross sections calculated from the next-to-leading order potentials 
is 8.6% at the peak and 2.4% at E — 5 GeV for the same value of a s (M|), the cross sections 
have come closer only slightly after the inclusion of the second-order correction to the potential. 
The remaining difference is much larger than what one would expect from an 0(a^) correction 
relative to the leading order, which is not fully included in our analyses, even if we take into 
account the high sensitivity to the coupling, cr pca k ex The purpose of this paper is to 

understand the origin of this unexpectedly large difference. 

As already mentioned, the difference of the cross sections reflects the difference of the poten- 
tials. The derivative of the potential dV(r)/dr is directly related to the size of the cross section; 
the cross section is larger if dV(r)/dr (= magnitude of the attractive force) is larger. This is 
because, with increasing probability that t and t stay close to each other, the wave function at 
the origin |^(0)| 2 increases, and so does the total cross section. Certainly, adding a constant 
to V(r) does not affect the size of the cross section at the peak. Thus, we Fourier transformed 
Vjkpt numerically from momentum space to coordinate space and plot the derivatives of the 
potentials in Fig. |2](a). To demonstrate the difference of the attractive forces, we show the 
difference of the derivatives of the two potentials, 

= — r dr ' (2) 

(solid line) in Fig. 0(b). We confirm that AF(r) > holds in the region probed by the 
toponium states, r ~ 0.03-0.1 GeV -1 . One also sees that both potentials have a common slope 
at r > 0.4 GeV -1 because of the severe constraints from the bottomonium and charmonium 
data. The kink seen in the figure is due to a discontinuity of g? 2 Vsfhmn/<^" 2 located at the 
continuation point, r = r c . 

In order to compare the asymptotic behavior of the potentials more clearly, we plot in 
Fig. ^|(a) the coordinate-space effective couplings defined by 

a v {l/r) = (-C F /r)- 1 V(r) (3) 

for Vjkpt ( r ) and VsfhmnM as solid and dashed lines, respectively. Contrary to our expectation, 
the difference of the couplings exceeds 3% even at very short distances, 1/r ~ 100 GeV. 

Naturally the question arises: Why is there such a large discrepancy between the potential 
constructed in momentum space and that constructed in coordinate space? To answer this 
question, let us examine a relation connecting the effective coupling in coordinate space, defined 
by Eq. (|3|), and the effective coupling in momentum space, defined from the momentum-space 
potential as 

a y (q) = (-4nC F /q 2 y 1 V(q). (4) 

The relation is derived from the renormalization group equation of oey(q) and exact to all 
orders. In the asymptotic region where the couplings are small, it can be given in the form of 
an asymptotic series frl3l, which reads numerically 



av(l/r) = «v + 1-225 a% + 5.596 c$ + 32.202 afy + ... (5) 



' In this paper we are not concerned with those differences of the cross sections which can be absorbed 
into an additive constant to the potential V(r), or equivalently, into a redefinition of the top quark mass. The 
theoretical uncertainty in the pole mass for our problem is discussed in 113] . 
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Figure 2: (a) Comparison of the derivatives of the potentials vs. r for (M§) = 0.118: 
dVjKPT /dr (solid line) and dVsFHMN /dr (dashed line), (b) Difference of the derivatives of the 
potentials vs. r. The solid line shows AF(r) = dV^w/dr — g?Vsfhmn / dr , and the dotted line 
shows AF(r) = dVj KPT /dr — dV ncw /dr. 
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Figure 3: (a) Comparison of the coordinate-space effective charges defined from Vjkpt (solid) 
and from Vsfhmn (dashed line), (b) Comparison of the coordinate-space effective charges 
defined from the various terms of Eq. (Eh. See the text for the description of each curve. 
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for rif = 5. On the right-hand-side, ay — oty{q — e~ 7B /r). All terms which are written 
explicitly are determined from the known coefficients of the f3 function, 0q , (5\ , and 0% . At 
present, we can use the above relation consistently only at 0(a\) because we know the effective 
couplings only up to the next-to-next-to-leading order corrections in perturbative QCD, i.e. we 
know the relation between ay and a^ only up to 0{a^). Due to this limitation, essentially, 
the effective coupling a v defined from Vsfhmn is the right-hand-side of the above equation 
truncated at the 0(a\) term, while ay defined from Vjkpt is the right-hand-side including 
all terms. Numerically, the 0{oty) term and the 0{a\) term contribute as +1.4% and +1.1% 
corrections, respectively, for ay = 0.1379 (corresponding to 1/r = 100 GeV). Therefore, these 
higher order terms indeed explain the difference of the effective couplings at small r. Fig. 0(b) 
shows several curves derived from the above relation: 

1. The solid line is ay {1/r) defined from Vjkpt- 

2. The dashed curve is ay + 1.225 a v , where ay = ay{q = e~ lE /r) is calculated using 
the perturbative prediction in momentum space. This curve is essentially the same as 
ay {1/r) defined from Vsfhmn ( r )> since it is the next-to-next-to- leading order perturbative 
prediction for the coordinate-space coupling at short distances. 

3. The dotted curve includes the next correction, 5.596 ay, which is in fact even larger than 
the 0{a\) term below 1/r ~ 30 GeV. 

4. The dash-dotted curve includes the 0{a Y ) term. 

We observe that the agreement of both sides of Eq. (|) becomes better as we include more 
terms at small r, while it becomes worse at large r on account of the asymptoticness of the 
series. From the purely perturbative point of view, the discrepancy between our two potentials, 
Vjkpt an d Vsfhmn, m the asymptotic region thus seems real, an indication of large higher order 
corrections. When the third-order correction to the potential will be computed in terms of a^rg 
in the future, the 0{ay) term will be treated consistently and the difference will reduce by 
1.4% at 1/r ~ 100 GeV. 

We may consider this difference of ay{l/r) at short-distances as an estimate of the higher 
order corrections on the basis of the following observations. First, the two potentials are equal 
up to the next-to-next-to-leading order, and there seems to be no reason a priori for considering 
one of the two to be more favorable theoretically. Secondly, if we apply the same method (the 
relation between dy and ay) to estimate the size of the already known O{o i s ) correction, we 
obtain 7t 2 /3q/3 = 193.4, which turns out to be a slight under-estimate of the true correction 
a 2 = 333.5 pj (n f = 5).| 

Moreover, the above 3% uncertainty of ay {1/r) at 1/r ~ 100 GeV provides a certain 
criterion for the present theoretical uncertainty of the ti cross section. In fact, it would already 
limit the theoretical accuracy of ay{l/r) at longer distances to be not better than 3%. If we 
combine this with a naive estimate cr pe ak oc a v , we expect a theoretical uncertainty of the peak 

J We may compare the coefficient of each color factor and find a certain similarity in them: 

7r 2 /?2/3 = 44. Ca - 32. C A T F n f + 5.8 {T F n f ) 2 

a 2 = 70. Ca 2 - 44. C A T F n f + 0.90 C F T F n f + 4.9 {T F n f ) 2 ' 
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cross section to be 5o" pea k/o" pea k ^ 6%. Therefore, the large discrepancy of the cross sections 
which we have seen turns out to be quite consistent with this estimated uncertainty. 

One is tempted to include one more term of the series @ to define a (new) coordinate-space 
potential despite our ignorance of the corresponding terms in the relation between ay and a^jg, 
since this would apparently reduce the difference between the two effective couplings. In fact 
we did this exercise, but (to our surprise) it did not bring the cross section closer to the one 
calculated from the momentum-space potential Vjkpt m the peak region. This cross section 
calculated from the potential V new (r), which incorporates the 0(a v ) term of Eq. (H), is shown 
dotted curve in Fig. 

We may understand the reason why the cross section did not approach that of Vjkpt if we 
look at the difference of the "forces", AF(r) = cIVjkpt / dr — dV new / dr , shown as a dotted line in 
Fig. ^|(b): it can be seen that, upon inclusion of the 0(a v ) term, the difference AF(r) decreased 
at small distances, r < 0.05 GeV -1 , as expected, whereas AF(r) increased at distances r > 
0.05 GeV -1 which is still in the range probed by the toponium states. It is due to a compensation 
between the decrease and increase of AF(r) that the normalization of the cross section scarcely 
changed. The increase of AF(r) at large distances results from the bad convergence of the 
asymptotic series, Eq. (||), for a large coupling, as we have already seen in Fig. |](b). This fact 
indicates that we are no longer able to improve the agreement of the cross sections by including 
even higher order terms, as we are confronting the problem of asymptoticness of the series. 

Some indications can be obtained by looking into the nature of the perturbative expansion 
of each potential. Within our present knowledge of the static QCD potential, the perturbative 
series looks more convergent for the momentum-space potential than for the coordinate-space 
potential. To see this, one may compare the (3 functions of the effective couplings (the V-scheme 
couplings) in both spaces |TD| . Numerically, the first three terms in the perturbative expansion 
read 

• (momentum-space coupling) 

^ 2 dav = _ QmQla 2 _ 0.2449 c4 - 1.198 <4 + ... (6) 

(JjLL 

• (coordinate-space coupling) 

^ 2 dav = _ Q 61Q1 -2 _ Q 2UQ -3 _ L945 ^ + . . . (7) 

(JjUj 

for rif = 5. The first two coefficients are universal. The third coefficient depends on the scheme 
(the definition) of the coupling. As the third coefficients for the V-scheme couplings are quite 
large, the third term of the (3 function is comparable to the second term already for ay = 0.20 
and for ay = 0.13, respectively!]] The difference of the third coefficients between momentum 

§ The shift of the peak position to lower energy is caused mostly by a decrease of the constant Co in Eq. ( |l3| ) 
and not due to an increase of the attractive force. Since the effective coupling ay(l/r) runs faster for Vnew, the 
perturbative potential is connected to the intermediate-distance phenomenological potential at a deeper point. 

» This is the reason why we evolve the MS coupling instead of evolving the V-scheme couplings using their 
own (3 functions. Otherwise we would have lost the reasoning to keep the third term of the (3 function at a 
fairly large momentum/short distance. For comparison, the (3 function of the MS coupling for the same n/ is 
given by 

dev. — 

-ffi = -0-6101 _ 0.2449 o^gg - 0.09116 + . . . . 
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Figure 4: Comparison of the cross sections for different choices of the scale. The upper 
three curves are for the momentum-space approach: /x = q (solid), /x = \[2q (dotted), and 
fj, = q/V% (dashed line). The lower three curve are for the coordinate-space approach: \i = 
fi2 = exp(— Je)/i" (solid), /x = \/2/X2 (dotted), and /i = (dashed line). 



space and coordinate space originates from the 7t 2 /3q/3 term in Eq. (0), which comes from the 
Fourier transformation. (Compare Eqs. (|T0| ) and (0).) Although the magnitude of the third 
coefficients is of the same order, in practice it makes a certain difference whether an apparent 
convergence is lost at ay = 0.20 or ay = 0.13 because there is a large scale difference between 
the two values. This indicates a worse convergence in coordinate space than in momentum 
space. 

If we evolve the coordinate-space coupling a v using its own (3 function up to the third term, 
the coupling exhibits an infrared pole at 1/r = A ~ 2.5 GeV, which is an order of magnitude 
larger than Aj^g of the MS coupling. The asymptoticness of the series in Eq. (H) is closely 
related to the existence of this pole. In fact, one may estimate the uncertainty caused by the 

Ar + (Ar) 2 



asymptoticness of the expansion to be 8a\{l/r) ~ Ar + (Ar) 2 + . . . fll3| . If we translate this 
to the uncertainty in the slope of the coordinate-space potential, we obtain 8F(r) ~ A 2 . This 
is in good agreement with the discrepancy AF(r) ~ 1-6 GeV 2 in the region r > 0.05 GeV -1 
in Fig. ||](b), where the usability of the asymptotic expansion is already limited to the first two 
or three terms. (For r < 0.05 GeV -1 , one may reduce the difference by including more terms.) 

It is interesting to examine the level of uncertainties within the momentum-space approach 
or the coordinate-space approach by itself. Fig. ^ shows how the cross section changes when 
we vary the scale by a factor of 2 in each approach: from [i — qf \/2 to /x — \/2q in Eq. (28) of 
PTOj] in the momentum-space approach (upper three curves), and from /x = /i 2 / y/2 to /x = y/2fi 2 
in Eq. (44) of [[L(J in the coordinate-space approach (lower three curves). For the momentum- 
space approach, the variation of the cross section is 2.2% at the peak and around 0.6% for 
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cm. energies above threshold. Meanwhile in the coordinate-space approach, the variation of 
the cross section amounts to 0.7% at the peak and 0.9% at larger cm. energies. These results 
may be regarded as an internal consistency check for each approach and even as a sign for the 
stability of the theoretical predictions. Nevertheless one should keep in mind that the internal 
consistency is not the same as the accuracy of the theoretical predictions. Since the 0(a 2 s ) 
corrections to the potential resulted in an unexpectedly large modification of the total cross 
section (see Figs. || and |6|), it would be legitimate to consider each of our results as accurate 
only if the same method could estimate the size of the next-to-next-to-leading order correction 
reasonably well and hence if the cross section became considerably less sensitive to the scale 
variation after including this correction. This is not the case in our problem, however. The very 
existence of a large constant at the next-to-next-to leading order (02 in |T0[ , see also Eqs. ( PH ) 
and (0)), which generates these large modifications, may indicate also large corrections at even 
higher orders. 

Still there may be some possibilities to reduce the difference between the momentum-space 
potential and the coordinate-space potential in the region probed by the toponium states. An 
obvious point to be improved is to remove the discontinuity of V' s / p HMN (r) at r = r c . If we 
employed a smoother interpolation of the perturbative potential to the intermediate-distance 
potential, we would have a better agreement of the two cross sections. This tendency is expected 
due to the specific interpolation method adopted for Vsfhmn(?")- It matches the perturbative 
potential exactly in the short-distance region up to a vicinity of the infrared pole. The rapid 
acceleration of the running of av(l/ r ) towards the pole tends to amplify the deviation from 
Vjkpt(^)- Therefore, if we employed a smoother interpolation to a phenomenological potential 
while keeping the potential to approximate the perturbative potential at short distances, the 
running of o:v(l/r) should be tamed, and hence the potential should come closer to Vjkpt(^), 
see Fig. 0(b). This tendency has been seen [15) at the next-to-leading order in the comparison 
of the cross sections calculated from Vsfhmn(?") and the Strassler-Peskin potential |]. 

One has to be careful with this argument, however. The slope of the potential in the 
intermediate-distance region is fixed by experimental data, which correspond to one fixed value 
of a^jg(M|) (= the true value in nature). We are interpolating the prediction of perturba- 
tive QCD, which obviously depends on our input value of a^jg(M|), to a phenomenological 
potential, which is independent of it. This means, we do expect a non-smooth transition for 
any value of a^g(M|) different from the true value. Moreover, if we want to extract the value 
of a^jg(M|) by comparing the theoretical predictions to the experimentally measured cross 
section, the sensitivity to a^g(M|) decreases if the predictions depend on the way we perform 
the interpolation. Ideally we would want to have an intermediate-distance potential as the 
prediction of QCD — necessarily non-perturbative — for a given input value of a^g(M|). 

So far we have examined the difference between the momentum-space potential and the 
coordinate-space potential in detail, and we have taken this difference as an estimate of the 
theoretical uncertainty of the QCD potential. This estimation may, however, be somewhat 
misleading. One might well argue that the difference is an artifact of our inadequate use of 
the perturbative expansion in describing the potentials related by Fourier transformation. For 
illustration, let us consider a hypothetical case where we know cty(q) exactly. In this case, if 
we calculate ay(l/r) via (numerical) Fourier transformation of the momentum-space potential, 
in principle we can calculate ay(l/r) to any desired accuracy by investing more time. On the 
other hand, if we calculate ay(l/r) using the series on the right-hand-side of Eq. (|5|) naively, 
there is a limitation in the achievable accuracy because the series is only asymptotic. Certainly 
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such a limited accuracy does not reflect any theoretical uncertainty of ay(l/r). This nature 
should not be confused with our claim. We claim that there is a limitation within perturbative 
QCD in relating ay(q) or ay(l/r) to a^g(Mf), and estimate theoretical uncertainties in this 
relation using the difference of the potentials in the two spaces. 

Another question that may be asked in connection with our estimate of the higher order 
corrections to the QCD potential is: Why should Fourier transformation know anything about 
the higher order corrections? We do not know the answer to this question, but we can say at 
least the following. It is not only the effect of "pure" Fourier transformation. The effect of 
lower-order QCD corrections are included in this estimation through (3q , @Y an d 0% > which 
determine each term of Eq. (|5]) . 

It would be important to understand the problem of the difference in the potentials also 
in momentum space, at least as much as we do in coordinate space presently. We have not 
done this analysis yet because of the difficulty in the numerical Fourier transformation of the 
potential from coordinate space to momentum space. 

Let us comment on the relation between this work and the work which has been com- 



pleted very recently. Ref. |12| presents a fixed-order calculation (without log resummations), 



which includes the full second-order corrections to the tt threshold cross sections. While we 
include corresponding corrections only to the static QCD potential in our analyses, we also em- 
ploy the renormalization-group improvement and thus resum large logarithms in the potential. 
In this sense the two works are complementary to each other. Both works give the common 
conclusion that the second-order corrections to the ti threshold cross section are large and may 
indicate a poor convergence of perturbative QCD, although they are based on qualitatively 
different arguments. The full set of the fixed-order 0(a 2 s ) corrections to the cross section near 
threshold (|T^] are larger in size and even more scale dependent than the corrections to the po- 
tential alone. The theoretical uncertainty may therefore be larger than indicated by the study 
in our paper. 



Summary 

• There is a difference between the potential constructed in momentum space and that 
constructed in coordinate space even at a fairly short- distance, 1/r ~ 100 GeV. The 
difference can be understood within the framework of perturbative QCD. We already 
know that there is a large correction at 0{a\) in the relation between the two potentials, 
although a consistent treatment is not possible until the full 0(af) corrections to the 
QCD potential are calculated. 

• The above difference at short distances would limit the theoretical accuracy of the QCD 
potential at longer distances and thus provides a criterion for our present theoretical 
uncertainty of the ti cross section, <5cr peak /<7 pea k ~ 6%. 

• In addition, it seems that we are confronting the problem of the asymptoticness of the 
perturbative series in the calculation of the ti cross section, as the top quarks do not 
probe a region which is sufficiently deep in the potential. We may not be able to improve 
our theoretical precision even if the higher order corrections are calculated in perturbative 
QCD. 
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• We may, however, discuss which of the two approaches gives a more favorable result 
theoretically. Up to the second-order corrections, the perturbative series looks better 
convergent for the momentum-space potential than for the coordinate-space potential. 
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A The momentum-space potential 

It seems appropriate to describe the potential Vjkpt used in the present analysis in some more 
detail. This potential is very similar to the potential Vjkt described in || and used in all later 
numerical studies within the momentum-space framework. It includes, however, the next-to- 



next-to-leading order terms from ||10|| . The momentum-space potential can be written as 

UjKPT(q) = V (q cut ) ■ (27r) 3 <5(q) - 4ttC f ftjKP j (q2) . ( 8 ) 

The effective coupling ckjkpt is defined to coincide with the two-loop perturbative prediction 
for large momenta, to be Richardson-like for small momenta, and to simply interpolate between 
these two shapes in some intermediate range: 

{«v, P crt(q 2 ), |q| > qi = 5 GeV 

«Rich(q 2 ), |q| < Q2 = 2 GeV (9) 

«Rich(q 2 ) + («v, P crt(q 2 ) - aRich(q 2 )), q2 < |q| < qi- 

The intermediate regime is only introduced to obtain a smoother transition between the small 
and large momentum parts, respectively. 

The first difference between the updated potential Vjkpt and the former version Vjkt is the 
fact that we are now able to use the full two- loop expression for the perturbative part, 



a Y ^ 2 )=^ m (c l 2 )(l + a 1 ^^ + a 2 {^^) ), (10) 



with the coefficients a\ and a 2 given in fH|. As the 6-quark threshold is neglected, rif = 5 is 



set throughout in the evolution of the MS-coupling, which can now consistently be performed 
at three-loop accuracy. 

A Richardson-like behavior for small momenta is chosen since the Richardson potential 



is known to describe the charmonium and bottomonium spectra fairly well. A pure Richardson 
form, however, would lead to severe numerical problems. Hence the ansatz has to be modified 
slightly by introducing two "subtraction terms" , 

/ 9n 4tt / 1 A R Ar \ 

«Rich q 2 = tw ^ rr: — - -f + * 2 ii 

A)w = 3 ) VM 1 + -p-) ^ q 2 + qcut/ 



11 



with Ar = 400 MeV. The first subtraction regulates the divergent behavior for |q| — > 0, 
the second subtraction is designed to reduce the modification introduced through the first to a 
minimum. Without the second additional term, the linear part of the position-space Richardson 
potential would be removed completely, whereas with it the first subtraction is cancelled for 
q 2 S> q 2 ut , and thus a big part of the confining potential is kept. It thus seems desirable to 
choose the parameter q cut small, but evidently it cannot be put to zero to really recover the 
pure Richardson potential. However, the linear part of the potential plays practically no role 
for the it-system as will be demonstrated below. The exact value of q cu t is therefore relatively 
unimportant and the adopted value q cut = 50 MeV results in both numerical efficiency and a 
fairly good accuracy of the predictions. 

The constant Vo(q cu t) in Eq. (|8|) is to some extent an arbitrary parameter. Different choices 
of Vo(q cu t) reflect the ambiguity in the definition of the pole masses for confined quarks. For 
Vjkpt (q) the choice 

47rC F Al . . 

A)(n f = 3) q cut 

is used. It leads to a Richardson- like potential that depends only weakly on the parameter q cut 
and coincides with the true Richardson potential in position space in the limit q cut — > 0. With 
this potential one obtains for the pole mass of the b quark m;, = 4.88 GeV. The choice of Vq is 
the second difference to the potential used in earlier works, where the constant Vq was fixed by 
the condition VjktO = 1 GeV" 1 ) = -1/4 GeV leading to m b = 4.7 GeV. 

In Fig. [| predictions corresponding to different momentum-space potentials are shown for 
the total it-production cross section as functions of energy E. The solid line corresponds to 
the potential Vjkpt of the present article. The dash-dotted line corresponds to the potential 
Vjkt using the one- loop formula for av, pe rt, the two-loop evolution for the MS-coupling and 
fixing Vq through V(r = 1 GeV" 1 ) = —1/4 GeV. There are two differences between Vjkpt and 
Vjkt- First, the inclusion of the two- loop correction to the perturbative potential increases the 
strength of the attractive interaction between i and i, and thus leads to an increase in the cross 
section. This is nicely demonstrated by the dashed curve, which corresponds to the inclusion 
of the two-loop potential and the same choice of V as in Vjkt- Second, the modified choice 
for Vq leads to a small shift of about 300 MeV in the energy scale, which is just the difference 
between the two Vq. The dotted curve has been included to demonstrate that the it system is 
quite insensitve to the long-range part of the potential: this curve corresponds to the choice 
q cut = oo, i.e. to completely removing the confining part from the potential Vjkpt and setting 
Vq = 0. 



B The coordinate-space potential 

The short-distance part of the coordinate-space potential is given by the next-to-next-to-leading 
order static QCD potential in position space flTOf , whereas its form in the intermediate- and 
long-distance region is determined phenomenologically. We thus have 

V (r) = I Tpert(^) at r<r c , 

sfhmni ) \ c + Ci log(r/r ) exp(-r/ri) + ar at r > r c . y ' 
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Figure 5: Comparison of the total cross section (normalized to R) for tt-production as a function 
of E = v/i - 2m t for m t = 175 GeV, a m = 0.118, r t = 1.421 GeV, calculated from the 
momentum-space potentials described in the text. The solid line corresponds to the potential 
Vjkpt used in the present analysis and the dotted line to Vjkpt with the confining part removed. 
The dash-dotted line corresponds to the potential Vjkt of ref.p and the dashed line to the 
inclusion of the two- loop correction to its perturbative part. The dashed line differs from the 
solid line by a constant shift in E. 



Here, 



V peit (r) 



-a 



1 + ai «Ms(^) 



4tt 



+ U 2 



)( 



MS 







47T 



(14) 



represents the coordinate-space potential in the second scheme, \ii = exp(— 7s)/r. The coeffi- 
cients a\ and (X2 are the same as in the momentum-space potential, and (3q = (IICa — 4Tprif)/3. 
See Ref. for details.Q 

The values of the phenomenological parameters Tq, ri, a and c\ are taken from Ref. and 
are tuned to reproduce bottomonium and charmonium data well: 



r 
r\ 
a 
ci 



0.2350 
3.745 
0.3565 
0.8789 



GeV" 1 
GeV" 1 
GeV" 2 
GeV 



(15) 



' We evolve the MS-coupling Q!Tjg(/i) by solving the three4oop renormalization group equation numerically 



for a given initial value at fi 
used in JhJ. 



Mz, whereas an approximate solution to the renormalization group equation is 
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We fix Co and r c by requiring that both the potential VsfhmnM and its first derivative are 
continuous at r = r c . For example, r c = 0.2526 GeV -1 and c = —1.972 GeV for a^g(M|) = 
0.118. 

This potential is an improved version of the potential proposed in || by including the next- 
to-next-to leading order terms to the short- distance QCD potential. We compare the cross 
sections calculated from the present version and from the old version in Fig. ||] for a^jg(M|) = 
0.118. 



R(E) 




E [GeV] 



Figure 6: Comparison of the total cross section (normalized to R) as a function of E = ^—2mt 
for m t = 175 GeV, c^jg = 0.118 and and T t = 1.421 GeV using the old and the present versions 
of the coordinate-space potential Vsfhmn- The dash-dotted line corresponds to the potential 
described in 0. The solid line shows the prediction of the present potential, Eq. (p~3|). 
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